Magnetic susceptibility of the quark condensate via holography 
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We discuss the holographic derivation of the magnetic susceptibility of the quark condensate. It is 
found that the susceptibility emerges upon the account of the Chern-Simons term in the holographic 
action. We demonstrate that Vainshtein's relation is not exact in the hard wall dual model but is 
fulfilled with high accuracy. Some comments concerning the spectral density of the Dirac operator 
are presented. 
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I. INTRODUCTION 

The analysis of the QCD properties by holographic 
methods is one of the most promising approaches to the 
description of the strong coupling region. The unique 
holographic model for QCD has not been found yet hence 
there is no hope to get the generic quantitative predic- 
tions at present. However there are some QCD results 
which seem to be independent on the details of the dual 
geometry hence one could consider these universal ob- 
jects or relations to test the holographic picture. On the 
other hand it is instructive to analyze if some relation 
is universal indeed testing it in the different holographic 
geometries. 

The simplest relation to be tested is the Gell-Mann- 
Rennes-Oaks one which was shown to be true in all holo- 
graphic models of QCD like hard wall models [ll, [20j . 
soft wall model [1] or Sakai-Sugimoto model Q- The 
main focus in our paper is the magnetic susceptibility of 
the quark condensate describing the response of the QCD 
vacuum on the external magnetic field. It was introduced 
in Q in the context of the sum rules and investigation 
of its numerical value was performed in [Ttt^[T7|. 
More recently using QPE arguments Vainshtein [7] ob- 
tained the expression for the susceptibility in terms of 
the known QCD quantities . However the status of this 
relation is questionable since both sides of the correspon- 
dence have different anomalous dimensions and it is not 
clear if the higher states could influence the answer, (see 

MM) 

In this paper we shall analyze the magnetic suscepti- 
bility in the holographic setting and shall focus mostly at 
the simplest hard wall model [Jllj (introduced in [2l| - [23| ). 
We shall consider the calculation of the three-point func- 
tion similar the two-point calculations in |^ and form- 
factor calculations in [HI ■ We shall consider the spe- 
cial kinematics of the three-point function related to the 
susceptibility. It turns out that the only nontrivial con- 
tribution to the correlator comes from the Chern-Simons 



term in the dual action and substituting the solutions to 
the classical equations of motion we get the result for the 
susceptibility which is close to the Vainshtein's relation. 

The paper is organized as follows. First in Section 
II we calculate the three-point function and the mag- 
netic susceptibility of quark condensate in the AdS/QCD 
hard wall model. In Section HI we survey some other 
approaches to the calculation of this value, namely via 
Chiral Perturbation theory and via the relation with the 
Dirac operator spectrum density. The conclusion is given 
in Section IV. To make the paper self consistent, we state 
in the Appendix some results of which we will use. 



II. HARD WALL ADS/QCD MODEL. 

A. Chern-Simons action and 3-point function 

Our aim is to calculate the correlation function of two 
vector and one axial currents. To make the holographic 
calculation we take the simple "hard wall" AdS/QCD 
model [ll,[8|. We will work in notation of @ and use some 
results, calculated in [l],!^ (see also Appendix). The holo- 
graphic action involves kinetic and Chern-Simons term : 

S = Sym{Al.Ar) + Scs{Al) - Scs{Ar) 



where 



ScsiA) - 



Nc 
247r2 



Tr I AF^ - Ia^F+—A'^ 
' 2 10 



It is clear that only terms containing 3 gauge fields 
{A, V, V)a.Te relevant for the calculation of correlator: 

Tr{ALFLFL - ArFrFr) ^ 

> 2Tr{VFvFA + VFaFv + AFyFv) 

AW 

The classical solutions for the fields have the form: 
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V;^{z,Q) = V,%Q)-V^,{Q\z) (1) 
Al{z,Q)^Al{Q)-A^,{Q\z) 
A^,{Q\ z) = P{Q)i,a^{Q\ z) + P{Q)l,a\\ {Q\z) 
V^,,{Q\z)^P{Q)l,,v{Q\z) 
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where V^^iQ) and Af,{Q) provide the sources for the oper- 
ators. We assume the currents to have arbitrary charges 
with respect to SU (2) group, so the gauge group struc- 
ture of the result is: 



(AVV) = 
Nc 



SASVSV 



Scs = 



127r2 



(TATvTy) ■ {VFvFa + VFaFv + AFyFy) 
-{TATvTy) ■ {VFvFa + VFaFv + AFyFy) 



We will work in the gauge A^ = Vz 
the integral over z, we rewrite it as: 



VFvFa + VFaFv + AFvFv ^ 



0. To figure out 



Ae^^'^p'^ J dzA^[d,V,Vp - V.d^Vp] 

+A^[dzV,d^Vp - dM.Vp] 
+d^Ap[V,dzVp ~ dzV^Vp] 



Substituting Fourier components of fields ([T]) and in- 
troducing the tensors 



Ptp{ki)P^Ak2)P^p{h) 



:ll 



we get (denote v — v{ky,), v = dzv) : 



{A±MVu{k2)Vp{k3)) = 
Nc 



X J dz (ik2cr)[aj_vv + a±vv ~ 2a±vv] 

— {ik3g.)[d±vv + a±vv — 2a±vv] 

We can add a surface term {—ik2 + ik3)dz{a±vv) in the 
action, in order to make the 3-point function vanish, if 
one of vector momenta tends to zero. This will lead us 
to the expression: 



Similarly 



(A±p{ki)V4k2)Vpik3)) = -^{TATvTy) 8^{ki + fc2 + fc3)e^.p. / dz {ik2„)avi - {ik3„)avv (2) 



{A\\p,iki)Vu{k2)Vp{k3)) = -^{TATvTy) 5'^{ki + ^2 + A;3)ell^p^ j dz {ik2a)avv - {ik3^)avv (3) 



r 



B. Solution for Aj_ 

Let us consider the equation of motion for A± at small 
similarly to calculation in The equation is Q: 



dz -dzAl + ^Al Al 







We will work with the bulk-to-boundary propagator 
aj^(z), which is defined in ([l} and denote R^g^h? = k^ = 
3 (|A.4[) . Using variable y — !j-z^ — a^z^, we get the 
equation 



4/3, 



+ 



2k^ma 

9a4 



9a2' 



-4/3, 



which is an inhomogeneous modified Bessel equation. We 
introduce here the dimension parameter a which equals 
^^^1/3 ^ 395^g„ (ggg (jX3l),(IX4|)- One can argue that 



the last term is negligible and solution to homogeneous 
part is 



(4) 



where constants are fixed by the conditions on the IR 
boundary y,n = a^z^i = 1.82 (see (|A.2[) ): 



dza{z)\. 



= 3ay^/'dya{y)\y^ 

^ = ^2/3(2/™); B = I_2/3{ym) 

and UV boundary: 







aiz)\z 



Fy^'^ By-i/3r(l/3) 



22/3 



F = 



22/3 



Br(i/3)' 



Given this solution (which corresponds to Q=0), we 
can compute /^r, using the recipe, described in Q (see 
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and q, assuming to be small enough. 



55 ^ 95 



(5) 



Solution for Ai 



To obtain the longitudinal part of the 3-point func- 
tion we need to find bulk-to-boundary propagator in the 
pseudoscalar sector. It is the solution to equations Q: 



' z 



(6) 



2,1 dy^/j 2 Pma 



-2/3 



9^ 



3 y 



^ + 0(9^g^), (8) 



where y = t^. 

The homogeneous solution, subject to the boundary 
conditions 7r(e) = 0, (/3(e) = l,9z7r(zm) = dz(p{zm) = is 
the same as for a± ^ as expected at = 0. The Green 
function of this equation is: 



G{u,v) 



yl/3^1/3 



x[Ch/3{v)+DK,/3{v)] 



(7) with C and D defined by the condition: 



where ip is related to the longitudinal part of as = 
d^ip. We introduce the function ijj{z) — ip{z) — 7r(z), and 
eliminate 7r(z) from the system to get an equation on 
with the dimensionless variable t = az 



It satisfies the equation 



tdt ( ia^v 



-dt^ = 0, 
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where q = /a^. Now we can substitute v{t) = -^t^ + 
—t and write down terms up to the first order in m/a 



1 



G{y,y')^S{y-y')-^. 



We can compute the correction due to the quark mass 
in ([8]). It is obtained by the integral: 



^(„)^^^^2fcW p y.l/3G(y,y'),y-2/3^(0)(^.) 



9a4 

2k^ma Fy^/^[Ah,^{y) + BK^^jy)] 
9a4 AD - BC 

p,,i/3r 



y'^'''[Ch/^{y') + DK„^[y')][Ah/:,{y')+BK„^{y')] 



9a4 



AD-BC 



Using the value of ym = 1.82, which corresponds to 
the IR boundary in our model, we get for the correction: 

V'^™^?/) = -1-0004 • 



9a2 



Consequently, the solution for ip with correction due 
to the quark mass is: 



^{z)^i:{z) + 0{q) 



= Faz[Ah,^{y) + BK^/^{y)]- 



9 a2 



-z^ (9) 



we neglect here the correction due to the q term assuming 
that Q < 3a = llSOMew. 



D. Magnetic susceptibility of the quark condensate 

In this Subsection we calculate the magnetic suscepti- 
bility X of the quark condensate defined as 



{q<7t,uq)F = x{QQ)Fi^f^ 



(10) 



In order to find magnetic susceptibility , we study the 
3-point function 

{Ai\^{-Q)V,iQ - ks)Vp{ks)) 

in the limit fca — ?> 0, according to |7:] where the following 
expression for the susceptibility has been obtained 



^ 8^2 J2- 



(11) 
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Consider the classical solutions for the vector fields, 
calculated in [H, Q . 

v{Q, z) = Qz (k,{Qz) + f°/^""^ /i(gz)) — ^ 1 

v{k, z) = k'^z (xaikz) + ^^^^I„{kz)) - 



and substitute them into the correlator ^ : 

= -^i^ATvTy) el^p^iiks^) j dz a\\{z)v{Q,z) 



TT 



^{TATvTy) el^p„{ik3„) 



X / dz L(o)(z) + <^("')(z)l [Q^zKoiQz)] 



In this integral due to the fast fall of the vector prop- 
agator Kq{Qz) we can take the boundary value of 
in the first term. The second term can be calculated 
explicitly 

dz^<-''\z)[Q^zKoiQz)] = (^"(0) [ dzQ^zKoiQz) = i 



dzip''"'\z)[Q^zKoiQz)] 



'dzz^lQ^zMQz)]^ 1.075^ 



9 a 



where we've used the result ^ and the relation (|A.3p . 
Finally, we get the expression for 3-point function with 
corrections: 



which can be matched the OPE of 

I 



27r2 



TVc Amf{qq)x 



0{1/Q') 



This comparison allows us to determine the magnetic sus- 
ceptibility of quark condensate x 



X = -2.15-^-2 



(12) 



in close agreement with the result of Vainshtein ([TT]) . 
This agreement is parametrical, but not numerical, be- 
cause due to the small value of /tt® in our model, we 
get Xmod = 11-5 Gev~^ which is too large. Anyway, tun- 
ing the parameters of the model (mainly (qq)) can help 
fix /tt to its real value and get the reasonable numerical 
agreement with Vainshtein's Xvain — 8.9 Gev~^. We've 
checked, that such variation of parameters do not affect 
coefficient in significantly. Namely, its value changes 
less than 5%, then the parameter ym which is propor- 
tional to (qq) /mi is varied in the wide range from 1 to 8. 



III. 



OTHER APPROACHES 



A. 3-point function in ChPT 

We can compute the same 3-point function in the Chi- 
ral Perturbation theory ^13|] and compare the result with 



AdS/QCD. Note that the chiral Lagrangian is derived in 
the Sakai-Sugimoto model hence the comment below 
can be considered as the justification of the Vainshtein 
[3] relation. To obtain the {AW) correlator, we con- 
sider the Wess-Zumino-Witten term in ChPT action and 
turning on axial and vector external currents. 



Z^PT 



d xL2 + Zwzw 



L2 = — {D,UD'^U^), 



with Dn = dull — irM + iUh 



Z[U,l,r]wzw = 



487r2 



d^: 



:eM-/3 {W{U, Z, r)^--"^ - W(l, I, r)"^"^) 
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[L ^ R). 



Here 



and {L ■(-> R) stands for 



The leading contribution is given by the tree diagram, 
including (ott) vertex from L2 and (ttvv) vertex from 
Lwzw 

AisM =V2F{df,na'') 



ALwzwIt^vv] 



'487r2 F 



X {2df,TTdyVaVp - df^TTVadyVp) 



We can check, that the longitudinal part is exactly the 
same as in [7], if we, formally, expand it in /Q^: 



{A\\,^{-Q)V,{Q)Vp{k)) 



Nc_ 
27r2 



7y2 



(Q^ is Euclidean momentum) 



B. Relation with the Dirac operator spectrum 

Let us also briefly comment on the different calculation 
of the magnetic susceptibility via the spectral density of 
the Dirac operator and introduce eigenfunctions of the 
Dirac operator in the external gluon field A 

D{A)ux{x) = \uy,{x) 

Then the standard definition of the spectral density reads 
as 



where V is Euclidean volume and the averaging over the 
gluon ensemble is assumed. The value of the spectral 
density at the origin is fixed by the Casher-Banks rela- 
tion while the linear term was determined comparing 
the different calculations of the correlator of the scalar 



currents [15|. In the perturbation theory the spectral 
density behaves as p{\) oc that is starting from the 
third order the universality is lost because of the mixing 
with the perturbative modes. We would like to note that 
the magnetic susceptibility is sensitive to the last "non- 
perturbative" quadratic term in the spectral density. 

To explain this point let us consider the " two-point 
loop diagram" with tensor and vector vertexes in terms 
of the eigenfunctions and eigenvalues of the Dirac opera- 
tor. The simple inspection shows that the susceptibility 
is expressed in terms of two different contributions. The 
first "diagonal" contribution reads as 



m dX 



P(A) 



(A2 + to2)2 



while the second " nondiagonal" contributions involves 
the following integrals 



d'^xux{x)x,^ux' {x) 



and double integrals over the eigenvalues J dX J dX' . The 
"diagonal" contribution is IR divergent and this diver- 
gence is expected to be canceled by the "nondiagonal" 
terms amounting to a kind of sum rules. On the other 
hand it is clear that quadratic term in the spectral density 
yields the finite contribution. It is not clear if the "non- 
diagonal" terms yield the IR finite contribution as well. 
This point does not allow us to write down the coefficient 
in front of the in the spectral density immediately. 
One can not also exclude that more careful treatment of 
the IR divergences should involve the derivation a kind of 
the effective action with the tensor insertion . We hope 
to discuss these issues elsewhere. 



IV. CONCLUSION 

In this paper we have derived the expression for the 
magnetic susceptibility of the quark condensate in the 
holographic QCD model. We have demonstrated that 
this object captures nontrivial anomalous properties of 
the dual model encoded in the Chern-Simons term. It 
vanishes if the CS term is not taken into account. The 
second important lesson concerns the validity of Vain- 
shtein's relation which is not exact but is fulfilled with 
the high accuracy. 

The numerical value of the susceptibility do not co- 
incide with recent estimations from the instanton liquid 
model [1^, [13] ,sum rules fit JJJ and phenomenology of 
D-meson decays [l9j. But it is calculated at significantly 
less energy scale: for our calculation Q <C 1 150 Mev, 
while others are calculated at Q ~ I Gev, so we do not 
find the contradiction. This also allows us to compare 
the result with Vainshtein, whose normalization point is 
about 0.5 Gev. 

The only parameter the coefficient in Vainshtein's re- 
lation depends on is the IR cut-off scale however the de- 



6 



pendence is very smooth. It would be interesting to dis- 
cuss the soft wall model and derive the dependence of 
the magnetic susceptibility on temperature and chemical 
potential. 



The boundary condition on (/) at z = e fixes the constant 
B': 

0(0, q) = Uq) = ~-^-Q2 B' + B' 
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Appendix 

In Appendix, we state some results of [J, 0] concerning 
the "hard wah" AdS/QCD model. The 5D coupling con- 
stant gs is fixed by the 2-point function of vector currents 
in til] ' 



i?2 



12^ 



(A.l) 



The position of the IR boundary Zm is related to the 
/9- meson mass 



B' 



-0o(g) 



therefore we finally get : 



dzcj){z) 



= g5 B Ml)— HQ e ) 

^WU=e =0 

We can compute the 2-point function of pseudoscalar 
currents, using the relation: 

which yields us the source for pseudoscalar current 

2mq(t) o (975(7) ■ 

In order to obtain the 2-point function, we vary the 
action twice with respect to 2mg0(O) and find: 



1 



Z2iMev 



(A.2) 



The parameter a is coupled with the value of quark con- 
densate (we take the value {qq) = (230Mev)^) and equals 



Nf{qq) 



{460Mev) 



(A.3) 



We shall also fix the constant A correcting calculation 
made in 0. First, compute the leading order solutions to 
the equation of motion for the pseudoscalar fields. These 
are solutions to the equations of motion ([5]), ([7]) with 
fixed boundary value oi (j) aX z = t. Differentiating ([7]) 
and substituting dz4> from ([6]) we get : 



v2\ z3 7;2 

^6 I 



Q'^dzTT 

z-^ 



dzTT = 0. 



We need to solve it near the boundary, so substitute 
asymptotic value v(z) = mz\z^o and denoting x — Qz it 
takes the form: 



0. 



At large we neglect the last term and obtain the mod- 
ified Bessel equation with A = . Hence the solution for 
7r(z) reads as : 

tt{z) ^ A'QzhiQz) -f B'QzKiiQz) - C . 
and using ([6|) we immediately obtain the solution for (j): 



R 

9l 



bd^,(j) — 



o*(«l+<?2)2 



JRQ 



2 r 



(2:, 91, 92) 



9l 



- A^R^m^ 



S(t){qi,z) 



dTT{qi,z) 



A^R^ 
dz(t>{q2,z) 



Sn—dzTr 
z-^ 



dzn{q2,z) 



hence, the pseudoscalar correlator is: 



{j:iq),Mq)') = 
= 26 



1 RQ' \ 55Wm2^,Q2 



9E 



-B'^^lniQ'e' 



m=Q 4 

Comparing with the QCD value [l^ : 



2 



{J^{q),Mqf) = S^'^QMQ'e') 



we find 



47r2 3g| 



fc2 = R^A^gl = 3 



(A.4) 
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